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I. INTRODUCTION 



Higher derivative operators produce negative and positive effects in quantum field the- 
ories. Among the negative effects are the lack of unitary of the S— matrix, the presence 
of negative-norm states (ghosts) and the violation of the Lorentz invariance. The ghosts, 
strictly speaking, become evident by reformulating a given higher derivative theory in terms 
of standard (lower derivative) operators, i.e., by removing from it the higher derivative op- 
erators by means of auxiliary fields. The lower-derivative theory obtained in this manner is 
a theory with indefinite metric where the "auxiliary" fields play the role of ghost fields. In 
addition, this lower-derivative reformulation of a higher derivative theory is a key step in 
its canonical quantization based on Ostrogradski's approach [lj. A positive effect, on the 
other hand, is the improvement of the ultraviolet (UV) behavior due to the exchange of both 
negative- and positive-norm states in the Feynman integrals. Another way to say this is that 
the propagators in theories with higher derivative kinetic operators are more convergent at 
the UV limit k 2 — >■ oo than in the usual ones, implying a better UV behavior of the Feynman 
integrals (see [2J, for a review on these issues). However, as reported in [3], such positive 
effect does not always occur due to subtle problems with the analytical continuation from 
the Minkowski to the Euclidean space. 

To circumvent or eliminate the negative effects and take advantage of the positive ones, 
it is important to recognize that they have its origins in the additional degrees of freedom 
introduced by the higher derivative operators. Therefore, to construct a quantum field 
theory with higher derivatives or with indefinite metric which satisfies the minimal physical 
requirements (unitarity, Lorentz invariance and positive-energy spectrum), it is necessary 
to devise suitable mechanisms to get rid of "runaway" solutions or troublesome degrees of 
freedom. Evidently, one means is to impose constraints or boundary conditions on certain 
sectors of the theory. 

A long time ago, Lee and Wick showed in |H |5] that it is possible to construct a quantum 
field theory with indefinite metric (or in another language, with higher derivative terms) 
in which the S— matrix is relativistic and unitary. Specifically, they proposed a variant of 
quantum electrodynamics (QED), which is the result of introducing "heavy" negative-norm 
fields (heavy ghosts) in the gauge and fermion sectors of the original QED theory, free of 
divergences. The unitary problem, as a result of the presence of ghosts in Lee- Wick QED, is 
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avoided by requiring (as a boundary condition) that ghosts do not belong to the asymptotic 
states of the S— matrix, a condition that is possible only if the ghosts are unstable particles 
(i.e., if they have a non- vanishing decay width). On the other hand, as pointed out by Lee 
and Wick [6j afterward, the violation of Lorentz invariance overlooked in Lee- Wick QED 
(see e.g. |2]) can also be avoided by adopting the prescription of Cutkosky et al. [7] in the 
choice of the Feynman contours. Nonetheless, in spite of these achievements, the Lee-Wick 
QED theory is still plagued by some residual ghost effects such as acausality [8]. 

Nowadays, higher derivative operators are common in several branches of quantum 
physics and their negative effects are adequatedly dealt with by applying the Lee- Wick's 
ideas or any other killing-ghost mechanism, for instance, the method of "perturbative con- 
straints" in nonlocal theories (see, e.g., the work by Simon in [9j). In lower-energy effective 
theories, for example, higher derivative terms occur naturally as a result of integrating out 
the massive fields of a more fundamental theory and truncating its perturbative expansion, 
while in gravity p]] such terms (quadratic or higher powers of the curvature tensor) are gen- 
erated dynamically by radiative corrections. Higher derivative operators were also studied 
in susy/sugra models, string theory, Randall- Sundrum models, cosmology, phase transitions 
and Higgs models, and in other contexts (see [21 QJ] and references therein). More recently, 
the Lee- Wick's ideas were applied in the framework of the standard model (SM) in order to 
solve the hierarchy problem. Indeed, it was shown in [12J that all quadratic divergent radia- 
tive corrections to the Higgs mass are completely removed by introducing higher derivative 
kinetic terms in each sector of the SM. 

In this paper we bring together the good features (in relation to the improved UV be- 
havior) of supersymmetry and higher derivative operators. In particular, we show that a 
single higher derivative kinetic operator inserted in the usual Wess-Zumino action is suf- 
ficient to remove, under a rather general assumption on the complex poles, all the susy 
remaining divergences of the two-loop scalar self-energy. Notice that in this kind of the- 
ory (where supersymmetry and higher derivative properties are combined) two completely 
different mechanisms of removing UV divergences are involved. The cancellation of UV 
divergences in higher derivative theories occurs due to the exchange of normal and ghost 
states, while in susy theories this cancellation is achieved by the exchange of virtual particles 
with opposite statistics. Here we show explicitly in connection with the two-loop self-energy 
how both mechanisms work together to give a finite result. 
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On the other hand, to obtain some insights into the vacuum structure of our higher 
derivative Wess-Zumino model, we compute the effective potential at one-loop in the super- 
field formalism. By implicitly imposing the susy condition of non-negative energy to throw 
away "runaway" solutions, we find that supersymmetry remains intact at one-loop order, 
while the rotational symmetry is (spontaneously) broken iff a specific condition for o\ ^ 
is satisfied. 

Our paper is organized as follows. In Sec. [IT] we discuss in general terms the three dimen- 
sional Wess-Zumino model and define our higher derivative Wess-Zumino model (HWZ 3 ). 
This model results of introducing three types of higher derivative operators in the usual 



Wess-Zumino action. In Sec. Ill we compute the effective potential at one-loop and study 



its minima. As mentioned above, it is shown that supersymmetry remains intact at this or- 



der, while the rotational symmetry defines two phases of the theory. In Sec. |IV| we analyze 
the UV behavior of the scalar self-energy up to two-loops. Here we explicitly show how the 
mutual cancellation of UV divergences happens. Finally our conclusions are given in Sec |V} 

II. HIGHER-DERIVATIVE WESS-ZUMINO MODEL 

Our starting point is the three-dimensional M = 1 Wess-Zumino model, 

S = J d 5 z^-^D a $D a $ + m$$-U^ , (1) 

where $ denotes a complex scalar superfield whose Taylor expansion is 

$(x,9) = <p(x) + 6 a ij a (x) -6 2 F(x). (2) 

Here (p(x) and F(x) represent bosonic fields and ip a {x) represents a fermionic field. Through- 
out the paper, we shall adopt the notation of jl3j . 

N — — 

The superpotential U in general involves terms of the form (D a ) D ($$) 2 , where D a = 
da + iO^dap denotes the susy covariant derivative, No and N$ are respectively the number of 
susy derivatives (D a ) and the number of complex scalar superfields in a typical interaction 
vertex. 

Note that the rotational symmetry &'(z) = e*°$ (z), where a is a constant phase, and 
the Lorentz symmetry of the Wess-Zumino model Q restrict the values of Np and ]V$. In 
fact, Lorentz invariance requires an even number of susy covariant derivatives, with their 
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spinor indices completely contracted, while rotational invariance requires an even number 

of $ and $ superfields. Consequently, and iV$ must be even numbers. In addition, iV$ 

has to be greater than zero: iV$ > 2. 

By imposing the power-counting renormalizability condition, it is possible to find the 

general form of U . To this end, we compute the superficial degree of divergence a; of a 

typical Feynman diagram. From the free superpropagator 

_ n2 _ 2 

<* ^ d ) $ ("*> °')) = ~i (* - *) > (3) 

which is obtained by inverting the kernel of the Wess-Zumino action (flj, and taking into 
account the Grassmann reduction procedure |13| . one may show that 

- = /V + 2-|-f, (4) 

where / denotes the index of divergence, 

1=^-2, (5) 

V is the number of vertices in a typical Feynman diagram, E is the number of external lines 
and Ud is the number of susy derivatives transferred to the external lines 

Since the index I depends only on the form of the interaction vertex, a simple analysis 
of Q shows that the Wess-Zumino theory is renormalizable when the condition / < is 
satisfied. This in turn means that U can be of three types: 

U oc ($$) 2 , U oc D a §D a §, U oc $$, (6) 

thereby the only genuine vertex is the first one, while the others (kinetic- and mass-like 
vertices) can be re-summed in the quadratic part of the Wess-Zumino action Q and then 
completely absorbed by suitably redefining the superfield $ and the parameters of the theory. 

Our goal now is to extend the Wess-Zumino model ([T]) by adding to it higher derivative 
operators both at the kinetic and interaction parts. These operators must fulfill all the 
symmetries of the conventional theory, that is to say, they must be susy, Lorentz, and 
rotational invariants. Furthermore, all must be Hermitian. There is indeed an infinity of 
such operators, but in this work we shall restrict ourselves to the following ones: 



$ 2 (D 2 $) 2 + $ 2 ( J D 2 $) 2 , (7) 
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where the two formers modify the kinetic part of the Wess-Zumino model, while the latter 
modifies its interaction part. 

Including these higher derivative operators in (1), with U = g($$) 2 , our Higher- 
derivative Wess-Zumino model in three dimensions (HWZ3) is described by 



S = J d 5 z^-^D a <$>D a $ + aD 2 <$>D 2 <$>- ^D a <$>DD a <S> +m<$>$ - g (<$><$>)' 



+A 



(8) 



It is very simple to check out that the mass dimensions of the coefficients in this action 
are [a] = -1, [b] = -2, [m] = 1, [g] = and [A] = -2. Thus we shall take a = 1/M 
and b = 1/M 2 , where M and M are in principle arbitrary and different mass parameters. 
However, as we are regarding the higher derivative operators as the residual (low-energy 
limit) effects of an underlying fundamental theory, it is important to keep in mind that 
a, b and A must be very small compared with the original parameters of the theory. In 
particular, M and M must be of the Planck mass order. It should also be noted that the 
A interaction is non-renormalizable within the usual Wess-Zumino theory, according to the 
previous discussion. 

In terms of the component fields, the HWZ 3 Lagrangian is given by 



f-' ^bos £fer ^inti 



(9) 



where 



C bos = (pDip + m (<pF + TpF) + FF + a ((pDF + IpDF) + b (^C 2 <p + FDF) 



'/er- 



(10) 
(11) 



and 



int 



-2g 



ip 2 i) 2 + lp 2 ^ 2 + (p(p (cpF + pF + 2ip?p) 



+ 2A 



IpF (FF + 2i)id^) 



+ip 2 F 2 + if 2 (FUtp - {&tff) 2 ) + h.c. 



(12) 



Here the contraction of spinor indices follows the north-west rule (\) and the square of a 
spinor includes a factor of 1/2 in its definition. So ip 2 = \ip a ip a and ipidip = ip'^idj^ipp, for 
instance. 
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By setting up the F, F equations of motion of the free (i.e. switch off the g and A 
couplings) part of 

(aD + m) (p + (&□ + 1) F = 

(aD + m) Tp + (&□ + 1) F = 0, (13) 

one can easily see that F and F become dynamical fields only when b is different from 
zero. In fact, by taking b — in (|9|, F and F play the role of auxiliary Lee-Wick fields by 
introducing the higher derivative operator for the scalar (p field only in the so-called on-shell 
case. 



III. THE EFFECTIVE POTENTIAL AT ONE-LOOP 



In what follows we compute the one-loop effective potential in order to study the vacuum 
effect of the higher derivative operators. As is well-known, the classical potential is given by 
the negative of the classical HWZ3 action ^ evaluated at constant fields. Hence, writing 
the fields p and F as 

<P=A^Wi + i*i] F = -^={a 2 + in 2 }, (14) 

where <7j and 7Tj are real constant fields, the classical potential V c i is given by 

Vd = -- {o\ + ul) - (a x a 2 + 7Ti7r 2 ) [m - g {a\ + ul) + A (a* + tt 2 2 )] . (15) 

Since the effective potential V possesses the rotational SO (2) symmetry, a\ + m\ = 
exp(iu) ((Tj +i7Ti), that inherits from the classical HWZ3 action, we shall take 7Ti = in 



(15) to simplify the analysis of the vacuum structure of the theory. Doing this and solving 
the Euler-Lagrange equation dV c i/da 2 = 0, one finds two solutions for o 2 . 

± -1 ± VI + 12AA 
" 2 = ' (16) 

where A = go\ (a\ — r), with r = m/g. As we shall see below r is an order-parameter 
related to the spontaneous breakdown of the rotational symmetry at the classical level. 



Eliminating the "auxiliary" field a 2 of the classical potential V c i by means of (16), one 
obtains 



V ± = - 

d 108AV 2 L 



1 - 18AA ± (1 + 12AA) 



(17) 



Figure 1: Typical curves of and V7 . In V£ graph, with A = 0, the solid line represents the 
unbroken phase (r = 0), while the dashed line the broken phase (r ^ 0) of the rotational symmetry. 
In graph, the solid line exhibits the behavior of an atypical potential curve. 

In order to analyze the physical implications of this expression, let us consider and 
separately. Taking the positive potential V£ and expanding it around A = 0, 

V+ = 9 -^l(a 2 1 -r) 2 [l- 2AA + 9A 2 A 2 + O (A 3 )] , (18) 

we find out that Vj^ is analytic at A = and positive definite (Vj' > 0), as required by 
supersymmetric grounds. Fig. [I] displays two characteristic curves of for r = and r / 0, 
setting A = 0. They show us that there are two phases associated to the rotational symmetry: 
an unbroken phase of the rotational symmetry corresponding to r = and an spontaneous 
broken phase corresponding to r ^ 0. Classically, note also that supersymmetry remains 
intact at both phases. An analysis for the A ^ case leads us to the same conclusions. 
Now taking the negative potential Vj and its expansion around A = 0, 

K7 = " " g2 ° l { z " [1 " 2AA + 9A 2 A 2 + O (A 3 )] , (19) 

we see that in this case Vj is singular at A = and negative definite without a lower-energy 
bound (see Fig. 1). So, from a physical point of view, is the only acceptable potential 
that is positive definite and bounded from below. 

Next we shall compute the one-loop contribution to the effective potential by employing 
the steepest-descent method [15J or, as it is also known, the Jackiw functional method [16J, 
implemented in the superspace. For simplicity and technical reasons, from now on we will 
set b = = A. The one-loop contribution is given by 

V 1 = -—In Det (A - B) - —In DetC, (20) 



where v = J d 3 



A(z, z') 
B{z, z') 
C (z, z') 



[m + aD-g (3a 2 + vr 2 ) + D 2 + 2g {?>a x a 2 + ^2) O 2 ] S 5 (z — z 
J dzidz 2 H (z, zi) C~ x (zi, z 2 ) H (z 2 , z) , 
[m + aD-g [a\ + 3vr 2 ) + D 2 + 2g {a x a 2 + 3^^) 6 2 ] 5 



(21) 



5 >z - A 



with 



H {z, z') = [-2ga^ x + 2g (a^ + a 2 vr 2 ) 9 2 ] 5 5 ( 



z — z 



(22) 



The superspace functional determinants which appear in V\ are evaluated using the (- func- 
tional method |T71 H5] ; taking advantage of the rotational SO (2) symmetry (i.e. setting 
7Tj = 0). The result can be cast in the form 



Vi 



d 3 k 



In 



■ 2 (fc 2 ) z + a(3)fc 2 + /3(3,3) 
a 2 (fc 2 ) 2 + a(3) fc 2 + /3(3,0) 



+ ln 



,2 /j,2 



fc 2 ) z + a (1) fc 2 + /3 (1, 1) 
a 2 (£; 2 ) 2 + a(l)£; 2 + /3(i,0) 



, (23) 



where, defining u> (n) = m — n ■ go\ 



a in) = 1 — 2au (n) and (3 (ni, n 2 ) = a; (ni) + n 2 ■ (2gai<7 2 ) 



(24) 



Note that setting a = this result reduces to the usual case where there are no ghosts in 
the theory, since the (k 2 ) 2 terms vanish. 

In the following let us focus on the case with a / 0. To solve the above integrals we 
express each numerator and denominator as the product of two binomial factors of the form 
(k 2 + M 2 ) and employ the formulas in [20J. In this way, we get 



V 1 = i M * ( 3 ' 3 ) 3 - M * ( 3 > °) 3 + M * C 1 ' X ) 3 - Mi C 1 ' °) 3 ] 



(25) 



* = + !- 



where 



M± (ni, n 2 ) = J 1 — 2au (rii) ± \[\ — Aaco (wi) — n 2 ■ [Sa 2 ga^a^). 
yj2a v 



(26) 



Adding (|15|), with 7Tj = = A, and (|25j), we get the one-loop effective potential: 

a x a 2 (m - ga\) ~ ^ J2 [ M * & ^ ~ Mi & °^ + M * C 1 ' ^ ~ Mi ( X > ^ 



a 2 



*=+,- 



(27) 



At first glance, this result seems fairly intricate, however, considering the a 2 — linear approx- 
imation for the one-loop contribution, we will get some light about its physical implications. 
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In three dimensions M = 1 susy, this approximation is enough to analyze the possibility of 
susy breaking [T7]. Thus, the one- loop effective potential is given by 



V 



(To 



eff 



a x a 2 [m - go\ + X (a, 01)] 



where 



X(a, ax) 



9_ 

4-7T 



M_ (1, 0) - M+ (1, 0) 3 [M_ (3, 0) - M+ (3, 0)] 



I - Lf/^CL) v 71 -4aa;(3) 

After eliminating the auxilary field o"2 by means of its equation of motion, one gets 

V e ff = \o\ [m - go\ + X (a, a x 



(28) 



(29) 



(30) 



Notice that V e ff > and its minima occur at o"i = or ui 7^ 0, provided that <n is the solution to 
the equation: 

go\ = m + X (a, 01) . (31) 

This result means that susy is unbroken at one-loop order, while the rotational symmetry 
is spontaneously broken at this order only if there is a non-trivial solution a% 7^ for the 



condition (31 ). 

On the other hand, the series expansion for V e ff around a = 

r 2 



V 



eff 



g 2 a\ gaflm-;: 



2vr 2 a 2 



na 



+ \a\ (m - ga\) 2 + O (a 10 ) 



(32) 



shows us that V e ff is singular at a = 0. This result is very important and shows that a higher 
derivative term with a small coefficient, no matter how small it is, cannot be treated as a 
perturbation in the lower-derivative theory. Indeed, this result is well known in non-susy 
theories with higher derivative operators studied long ago in the literature |9]. 



IV. ULTRAVIOLET ANALYSIS OF THE SELF-ENERGY UP TO TWO-LOOPS 

In this section, we analyze the UV behavior of the scalar self-energy for the HWZ3 
model. We show that the introduction of a higher derivative operator in the Lagrangian 
(and consequently in the propagator) improves the structure of the divergences that appear 
in the momentum-space integrals, in contrast with the usual Wess-Zumino model. In fact, 
we shall see that all the integrals appearing in these corrections are finite in the UV regime, 
and so we have a much better UV behavior when comparing with the usual case. 
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Figure 2: Self-energy radiative corrections up to two-loops in the HWZ3 model. 

The Feynman diagrams which contribute to the scalar self-energy up to two-loops are 
depicted in Fig. [2] As in the previous discussion about the effective potential at the one- 
loop order, here we set b = = A as well. Hence the free propagator of the HWZ 3 model is 
given by 

A(k, 6 - 6') = (0\T${k,6)${-k,6')\0) 



1 



-\m-ak 2 -D 2 }5 2 (9-6') 



{m — ak 2 ) 2 + k 2 
t[A(k) + B(k)D 2 }5 2 (6-6') , (33) 



with 



A W = ( m _ a p)2 + k 2 > B \ k ) = ~( m _ ak 2y + k 2 ■ ( 34 ) 

In three dimensions and within the dimensional reduction (DReD 3 ) scheme [21 J , the 
one-loop diagram of Fig. [2] is finite. This statement is very general and well-known in lower- 
derivative theories [22] and is not spoiled by the introduction of higher derivative operators 
[14| . Indeed, as far as the UV behavior of higher-derivative theories is concerned, the Wick 
rotation can be performed in the usual way disregarding all contributions of the complex 
poles in the complex energy plane. The reason is that in the majority of the theories with 
higher derivative operators the residue contributions of the complex poles in the Cauchy's 
residue theorem are finite and so the form of the energy contour in the Feynman integrals 
is irrelevant. Note however that in order to construct a relativistic and unitary S— matrix 
it is necessary to fix the Feynman contour acoording to the Cutkosky's prescription [7]. 

Bearing these remarks in mind we proceed with the UV analysis of the diagram (b) in 
Fig. |2j The expression for this diagram is written as 

G 2b = c 2b J dpdkdq j d 2 M 2 2 $(0i,2#(0i, -p)A(fc, e 1 - e 2 )A(k, e 2 - ejAfa, e 2 - e 2 ) , (35) 

where p (with dp = d 3 p/(2n) 3 and so on) is the external momentum, k and q are the 
momenta appearing in the loops and C 2 b is a numerical constant and since we are concerned 
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only with the divergent behavior of the diagram we do not worry about it. After carrying 
out the .D-algebra, the expression above takes the form 



Q 2h oc J dpdkdq J d 2 9 {A{k)B{k)B{q)§{6, -p)$(9,p) + [B(k)] 2 B(q) [D 2 $(6, -p)Q(9,p) 

+D Q $(0, -p)D a <S>(6 lP ) + $(0, -p)D 2 <f>(e,p)] } , (36) 

and as can be easily seen the integrals in the momenta k and q are independent, and so by 
the same argument given above each one of these one-loop integrals is finite in the DReD3 
scheme. 

Now let us consider the last diagram in Fig. [2] Its mathematical expression is given by 



/2c 



c 2c dpdkdq d%d%^(9 1 ,p)^(e 2 ,-p)A(k,e 1 - e 2 )A(r,e 1 - e 2 )A( q , e 2 -e x ) , (37) 



where r = p + k — q and C 2c is a numerical constant and we do not worry about it, in the 
same way that we did with the constant C 2 b in diagram (b). After the .D-algebra, one obtains 



g 2c oc / dpdkdq J d 2 6 {[A(r)B(k)B(q) + A(q)B(k)B(r) + A(k)B(q)B(r)}<$>(9, -p)$(0,p) 

+B(k)B(q)B(r)D 2 §(6,-p)<$>(6,p)} . (38) 

The analysis of the UV behavior for this diagram is much more cumbersome. It presents the 
so-called "overlapping momenta", which did not appear in diagram (b). In order to isolate 
the divergent part from the finite part of the integral above, we must take its Taylor series 
expansion around p 2 = 0. Since each differentiation with respect to the external momentum 
improves the convergence of the integrand, the divergences will reside only in the first few 
terms of its Taylor expansion. In our case the Taylor expansion of (38) around p 2 = is 
given by 

02c oc JdpJd 2 6{n(p 2 )D 2 ^6,-p)^6,p)+S(p 2 )$(d,-p)$(6,p)] , (39) 

where 

K(p 2 )=n 1 +p 2 K 2 + (D(p 4 ), (40) 
S{p 2 ) = S l +p 2 S 2 + 0{p i ) , (41) 
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with, defining I = k — q and f(x) = (m — ax) 2 + x, 

Kl = S dkdd immiwy 

v f,r,~ r 2a l 2 + r 2b f(l 2 )+r 2c l 2 f(l 2 ) 

./ 1 f(k 2 )f(q 2 )[f(l 2 )} 3 ' 

Tj ~sia + s lb (k 2 + q 2 + l 2 ) 



Si = dkdq 



f(k 2 )f(q 2 )f(l 2 ) 



& = / dkdq [g2a(/2)2 + ( ' 26 + S2cP)ik2 + q2)]fil2) + [P(/2) + * ' gg( ^ 2 



/(A: 2 )/(g 2 )[/(/ 2 )] 3 

Here r« and s< denote numerical constants which might depend on the parameters a and m, 
while P and Q are polynomials in Z 2 of degrees 4 and 2, respectively. 

At this point two comments are in order. Firstly, since the two-loop self-energy in the 
usual Wess-Zumino model, i.e., without higher derivative operators, involves merely loga- 
rithmic and linear divergences (see Eq. Q, setting / = 0), we focus our attention only to 
the first term in the Taylor expansions for 1Z and S. Clearly, these terms enclose all the di- 
vergences present in the two-loop self-energy of the HWZ3 model. Secondly, we should point 
out that the (k, q)— integrals associated with p^— linear terms in the Taylor expansions are 
identically zero, a result that is in agreement with the Lorentz invariance of the model. 

As a result of the indefinite metric in Lee- Wick theories, in particular our HWZ 3 model, 
the IZi and Si Feynman integrals turn out to be finite. This assertion can be proved explicitly 
by evaluating each integral and observing the mutual cancellation between the divergent 
contributions from positive- and negative-norm (ghost) states, or implicitly by examining 
the superficial degree of divergence of each diagram. 

As the algebraic manipulations required to solve the IZi and Si integrals are very lengthy 
to be shown here, we are going to illustrate the divergence cancellation by evaluating in 
detail only the simplest integral, namely 1Z\. Up to a constant, we first rewrite the integral 
as 

f 7 1 

TZi = dkdq 



[(m - ak 2 ) 2 + k 2 }[(m - aq 2 ) 2 + q 2 }[(m - al 2 ) 2 + I 2 } 
d ~ kd \ 6 {k 2 + a){k 2 + f3){q 2 + a){q 2 + f3){l 2 + a){l 2 + J) ' (42) 



with 



a, P = -—^ (l — 2am =F \/\ — 4am) 
13 



Next we use the method of partial fraction decomposition to write 

1 1/1 1 



(k 2 + a){k 2 + /3) (3-a\k 2 + a k 2 + 
and split the integrand into eight fractions, obtaining 



Hi 



i 



a G ((3-a) 3 



dkdq 
1 



1 



(k 2 + a)(q 2 + a)(l 2 + a) (k 2 +a)(q 2 +a)(l 2 +P) 
1 1 



+ 



+ 



(k 2 +a)(q 2 + P)(l 2 + a) (k 2 + a)(q 2 + /3)(/ 2 + 0) (k 2 + f3)(q 2 + a)(l 2 + a) 
1 1 1 



+ 



(43) 



(k 2 +(3)(q 2 +a)(l 2 +(3) (k 2 +f3)(q 2 + f3)(l 2 + a) (F + /3)(g 2 + /3)(/ 2 + /3)_ ' 

Whitin the dimensional regularization scheme, each one of these integrals is calculated using 
the following formula [20J : 

d D kd D q 1 



X(x,y,z) 



(2n) 2D (k 2 + x)(q 2 + y)(l 2 + z) 



-2c 



327T 2 



7s + 1 - In 



ATTfi 2 



+ B r 



(44) 



where e = 3 — D, with D labeling the dimension of the spacetime, 7^ = 0.5772 is Eu- 
ler's constant and B res denotes the sum of residues of the integrand over all complex poles 
inside an energy contour appropriate for performing the "Wick rotation" (i.e. an energy 
contour which permits to change the real integration axis to the imaginary axis by means 
of the Cauchy's residue theorem). Strictly speaking, the analytical continuation from the 
Minkowski to the Euclidean space, which is common in conventional quantum field theories, 
is lost in the theories with higher derivative operators (i.e. Lee-Wick theories) due to the 
presence of complex poles in the complex energy plane. This fact is reflected, for example, 
in the non-vanishing value of B res that one finds in this sort of theory. 



Using the formula (|44|), we express 1Z\ in the form 

[X(a, a, a) - 3X(a, a, (3) + 3Z(a, /?, (3) - X(/3, /3, (3)} 



1 



a 6 (/3-a) 3 



(45) 



From (44) and (45), and assuming without proof that B res is finite, one can see that the 



divergent parts of the X terms cancel mutually, so that the integral 1Z\ as a whole is in 
fact finite, as was claimed before. The same procedure can be used to calculate the other 
integrals and all of them happen to be finite. The finiteness of the self-energy in the HWZ3 
model is in contrast with the usual Wess-Zumino model, in which the diagram (c) in Fig. 
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[2] gives a non- vanishing divergent contribution in the UV regime, showing us that a higher 
derivative kinetic operator improves the behavior of the model in this regime. 

There is a more elegant and general form to see why this better UV behavior is achieved. 



This is the study of the superficial degree of divergence w of a diagram. Recalling Eq. (33), 
we can see that the two terms in the expression for the superpropagator, one proportional 
to A(k) and the other proportional to B(k), give different contributions to a given diagram, 
since A(k) has power —2 and B(k) has power —4 in the momentum k. Moreover, B(k) 
is accompanied by a superspace derivative D 2 , which is not the case for A(k). We have 
therefore to consider that the complete propagator comprises two types of propagators in 
order to compute the superficial degree of divergence. These two propagators are defined by 



the A(k)— and B(k) — terms in (33): 



A(fc, 9-9') = A A (k, 9-9') + A B (k, 9-9'), (46) 

where 

A A (k,9-9') = i{m ~fl h2 s\e-e') 

[m — ak l Y + k z 

A B (k,9-9') = % °* Me-ff). 

(m — ak 1 ) 1 + k 

Taking into account the two types of propagators, it is straightforward to show that 

u = 2-x, (47) 

where 

x = 2V + P B + T f, (48) 

and, as before, V and rip are respectively the number of vertices and the number of susy 
derivatives transferred to the external lines, and P B is the number of propagators of the type 
B in a given diagram. Since x is strictly an integral positive number (x > 1), our HWZ3 
model exhibits only logarithmic (u = 0) and linear (00 = 1) divergences. 



On the other hand, from Eqs. (|36[) and (38) we see that rip assumes only two values 



rt£, = and rip = 2 and, for the two- loop diagrams in Fig. [2j we have that V = 2. Thus, 
by considering the worst case (i.e. P B = 0) in both diagrams, one sees that u> < 0, showing 
that these diagrams are finite. 
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V. SUMMARY AND CONCLUSIONS 



In this work we investigated the classical (and quantum) effects of three types of higher 
derivative operators introduced in the Lagrangian of the three-dimensional Wess-Zumino 
model. These operators respect all the symmetries of the original model, but the potential 
A— operator turns out to be non-renormalizable by power- counting arguments. At the clas- 
sical level, we show that these Lee-Wick operators modify the structure of the equations of 
motion for the component fields. In particular, one finds that the Lee- Wick b— operator 
promotes the component field F (P) from an auxiliary field to a dynamical field. 

We also considered the quantum aspects of the model in two distinct analysis. First, we 
computed the one-loop correction to the classical potential, which is singular at a null value 
for the higher derivative parameter a. This fact was already expected [9J and shows that the 
higher derivative term cannot be treated as a perturbation in the lower derivative theory. 
After this, we analyze the UV behavior of the self-energy up to two-loops. We showed 
by direct computation of one of the momentum-space integrals and also by calculating the 
superficial degree of divergence that the "setting-sun" diagram gives a finite contribution to 
the self-energy, in contrast with the usual Wess-Zumino model. This explicitly shows that 
the introduction of the higher derivative operator improves the behavior of the theory in the 
UV regime. 

As future efforts, we shall intend to analyze the features of the model presented here 
with b, A 7^ 0. This is a more involved work, due to the fact that the integrals appearing 
in the one-loop correction to the effective potential and in the correction to the self-energy 
are much more cumbersome and, besides that, the number of two-loop diagrams increases 
considerably. Also, we will consider the study of three-dimensional supersymmetric models 
with gauge fields within the framework of higher derivative models, a very interesting task 
which deserves attention. The situation in this case is different from what we have in the 
present work, for the gauge potential multiplet is a spinorial function on superspace and 
the procedure to obtain the propagators, even without higher derivative operators, is more 
intricated [18] . 
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